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ABSTRACT 

Recent observations suggest that neutron stars in low-mass X-ray binaries rotate within a narrow range of 
spin frequencies clustered around 300 Hz. A proposed explanation for this remarkable fact is that gravitational 
radiation from a steady-state r-mode oscillation in the neutron star's core halts the spin-up due to accretion. For 
the neutron star transients, balancing the time-averaged accretion torque with gravitational wave emission from 
steady-state, constant amplitude r-mode pulsations implies a quiescent luminosity too bright to be consistent with 
observations (in particular of Aql X-l). The viscous dissipation (roughly lOMeV per accreted nucleon for a spin 
of 300 Hz) from such an r-mode makes the core sufficiently hot to power a thermal luminosity ~ 10 34 ergs _1 when 
accretion halts. This is the minimum quiescent luminosity that the neutron star must emit when viscous heating 
in the core is balanced by radiative cooling from the surface, as is the case when the core of the star is superfluid. 
We therefore conclude that either the accretion torque is much less than M(GMR) X I 2 , or that a steady-state r-mode 
does not limit the spin rate of the neutron star transients. Future observations with Chandra and XMM promise to 
further constrain the amount of viscous dissipation in the neutron star core. 

Subject headings: accretion, accretion disks — stars: individual(Aquila X-l) — stars: neutron — stars: 
oscillations — stars: rotation — X-rays: stars — gravitation — gravitational waves 



1. INTRODUCTION 

With the launch of RXTE, precision timing of accreting neu- 
tron stars has opened new threads of inquiry into the behavior 
and lives of these objects. The neutron stars in low-mass X- 
ray binaries (LMXBs) have long been thought to be the pro- 
genitors of millisecond pulsars (see Bhattacharya 1995 for a 
review), and a long-standing observational goal has been the 
detection of a spin period of a neutron star in an LMXB. Re- 
cent observations (see van der Klis 1999 for a review) have fi- 
nally provided conclusive evidence of millisecond spin periods 
of neutron stars in about one-third of known Galactic LMXBs. 
Altogether, there are seven neutron stars in LMXBs with spin 
periods firmly established by either pulsations in the persistent 
emission (in the millisecond X-ray pulsar SAX J1808. 4-3658; 
Wijnands & van der Klis 1998) or oscillations during type I X- 
ray bursts (so-called burst QPOs, first discovered in 4U 1728— 
34; Strohmayer et al. 1996). There are an additional thir- 
teen sources with twin kHz QPOs for which the neutron star's 
spin may be approximately equal to the frequency difference 
(van der Klis 1999). A striking feature of all these neutron 
stars is that their spin frequencies lie within a narrow range, 
260 Hz < j/ sp i n < 589 Hz. The frequency range might be even 
narrower if the burst QPOs seen in KS 1731-260, MXB 1743- 
29, and Aql X-l are at the first harmonic of the spin frequency, 
as is the case with the 581 Hz burst oscillations in 4U 1636-536 
(Miller 1999). If this is the case, then the range of observed 
frequencies is 260 Hz < is sp i n < 401 Hz. The neutron stars in 
LMXBs accrete at diverse rates, from 10~ 11 MQyr" 1 to the Ed- 
dington limit, 10~ 8 MQyr _1 . Since disk accretion exerts a sub- 
stantial torque on the neutron star and these systems are very 
old (van Paradijs & White 1995), it is remarkable that these 
neutron stars' spins are so tightly correlated, and that none of 
the neutron stars are rotating anywhere near the breakup fre- 



quency of roughly 1 kHz. 

Observations therefore suggest that neutron stars in LMXBs 
are somehow stuck within a narrow band of spin frequencies 
well below breakup. Two explanations for this convergence of 
spin frequencies have been proffered. White & Zhang (1997) 
argued that the magnetospheric spin equilibrium model (see 
Ghosh & Lamb 1979 and references therein), which is appli- 
cable to the accreting X-ray pulsars, is also at work in LMXBs. 
In this scenario, the neutron star's magnetic field {B ~ 10 9 G) 
dominates accretion near the stellar surface, and the Keplerian 
period at the magnetospheric radius roughly equals the spin pe- 
riod, so that the accretion stream exerts no net torque on the 
star. Because the sources' luminosities (and presumably accre- 
tion rates) vary by several orders of magnitude, White & Zhang 
(1997) noted that this explanation requires either that the accre- 
tion rate be tightly correlated with the neutron star's magnetic 
field, B oc M 1 ' 2 , or that the torque be roughly independent of 
accretion rate when the magnetospheric radius approaches the 
radius of the neutron star. Moreover, the persistent pulses typi- 
cal of magnetic accretors must also be hidden most of the time. 

The other class of theories, first considered by Papaloizou 
& Pringle (1978) and Wagoner (1984), invoke the emission of 
gravitational radiation to balance the torque supplied by ac- 
cretion. Bildsten (1998) proposed that equilibrium between 
the accretion torque and gravitational radiation can explain the 
narrow range of observed spin frequencies. The source for 
the gravitational radiation could be a mass quadrupole formed 
by misaligned electron capture layers in the neutron star's 
crust (Bildsten 1998). Alternatively, as proposed independently 
by Bildsten (1998) and Andersson, Kokkotas, & Stergioulas 
(1999), current quadrupole radiation from an unstable r-mode 
oscillation (Andersson 1998; Friedman & Morsink 1998) in the 
liquid core of the neutron star could also limit the spin, as might 
occur in hot, newly born neutron stars (Lindblom, Owen, & 
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Morsink 1998; Owen et al. 1998; Andersson, Kokkotas, & 
Schutz 1999). Because the accretion rate of LMXBs does vary 
by several orders of magnitude, the small range of v spm among 
these objects also requires a correlation between the quadrupole 
moment and accretion rate. This correlation is much less re- 
strictive, however, than for magnetic equilibrium theories be- 
cause of the steep dependence of gravitational wave torque on 
the spin frequency. 

These theories have renewed interest in accreting neutron 
stars as gravitational wave sources. If gravitational radiation 
does in fact halt the spin-up of accreting neutron stars, then, 
regardless of the mechanism producing the gravitational radi- 
ation, the brightest LMXBs (such as Sco X-l, with dimen- 
sionless strain h c ^>2x 10~ 26 ; Bildsten 1998) are also promis- 
ing sources for ground-based gravitational wave interferome- 
ters, such as LIGO, VIRGO, GEO, and TAMA (Bildsten 1998; 
Brady & Creighton 1999). It is not certain, however, that ac- 
creting neutron stars in LMXBs do emit gravitational radiation. 
The only evidence to date is their narrow range of spin frequen- 
cies. It is therefore important to look for astronomical observa- 
tions, doable today, that can either corroborate or rule out the 
various mechanisms for gravitational radiation from LMXBs. 

In this paper we present a new observational test for r-mode 
driven gravitational radiation from neutron stars in one set of 
LMXBs, the soft X-ray transients. These are LMXBs in which 
accretion outbursts, lasting for days to months, are followed by 
periods of quiescence, lasting on the order of years to decades. 
Typical time-averaged (over the recurrence interval, rather than 
just over the outburst) accretion rates (M) for these sources are 
<; 10~ 10 Mq yr" 1 , smaller than those in the brighter persistently 
accreting LMXBs. We show that the quiescent X-ray lumi- 
nosities of these neutron star transients (in particular Aql X- 
1, which exhibits burst QPOs with a frequency 549 Hz; Zhang 
et al. 1998) can be used to determine whether r- modes with am- 
plitudes sufficient to balance the accretion torque are present in 
their cores. 

Recent theoretical (Brown, Bildsten, & Rutledge 1998) and 
observational (Rutledge et al. 1999b) works suggest that at least 
some fraction of the quiescent luminosity of a neutron star tran- 
sient is thermal emission from the neutron star's surface. Mo- 
tivated by the possibility of indirectly measuring the core tem- 
perature of an accreting neutron star, we consider the amount of 
heat that must be lost, on average, by the neutron star to main- 
tain a thermal steady state. If the spins of neutron star tran- 
sients are set by the equilibrium between the time-averaged ac- 
cretion torque and gravitational wave emission by steady-state 
(i.e., constant amplitude) r-mode pulsations in their cores, the 
required amplitude of the pulsations can be computed (§ fy. 
The steady-state assumption implies a certain magnitude of vis- 
cous dissipation, i.e., heat deposited directly into the core of 
the neutron star. If the core is superfluid, Urea neutrino emis- 
sion is suppressed and this heat escapes as thermal radiation 
from the surface of the star. We show (§ |]) that in this case 
the X-ray luminosity in quiescence, L q , would be about 5-10 
times greater than that observed. If the nucleons in the core 
are normal, then, as shown by Levin (1999), r-mode pulsations 
are thermally unstable (at least for saturation amplitudes of or- 
der unity). In this case it is unlikely that r-modes are currently 
excited in any of the known Galactic LMXBs. If for some rea- 
son a thermal steady state could be achieved in a normal fluid 
core, however, then Urea neutrino emission would carry away 
most of the r-mode heating, and the resulting lower quiescent 
thermal luminosities would be consistent, within uncertainties, 



with observations. Our test does not depend on how the r-mode 
is damped, but only on the assumptions that the dissipated en- 
ergy is deposited into the thermal bath of the star and that the 
star has reached a rotational and thermal steady state. We are 
only inquiring into total energetics, i.e., whether the viscous 
heating present matches that required by the spin equilibrium 
with the accretion torque. 

2. R-MODE VISCOUS HEATING OF ACCRETING NEUTRON STARS 

Recently Andersson (1998) and Friedman & Morsink (1998) 
showed that gravitational radiation excites the r-modes (large 
scale toroidal fluid oscillations similar to geophysical Rossby 
waves) of rotating, inviscid stars. Lindblom et al. (1998) com- 
pared the gravitational wave growth timescale r gr for the r- 
modes with the viscous damping timescale r v set by shear and 
bulk viscosities for normal fluids (i.e., no superfluidity); at ro- 
tation rates Vl ^ 0.065^, where tt K = (GM/fl 3 ) 1 / 2 is the Kep- 
lerian angular velocity at the surface of the star, the damping is 
sufficient to preclude unstable growth. The modes are excited, 
however, over a wide range of spin frequencies and tempera- 
tures that includes typical values for the neutron star transients. 

Gravitational waves radiate away angular momentum at a 
rate 
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is the canonical angular momentum of the (l = 2,m = 2) r-mode 
(Friedman & Schutz 1978; Owen et al. 1998), a is the dimen- 
sionless amplitude of the mode, and J is a dimensionless con- 
stant that accounts for the distribution of mass in the star (Owen 
et al. 1998). The gravitational wave growth time r gr is nega- 
tive, which implies instability. In a rotational steady state, this 
angular momentum loss is balanced by the accretion torque, 
A^ accr . For a fiducial torque, we assume that each accreted par- 
ticle transfers its Keplerian angular momentum to the neutron 
star, with a net accretion torque A^ accl - = M{GMR) X / 2 . Using r gl - 
as evaluated by Lindblom et al. (1998), we find the steady-state 
r-mode amplitude (Bildsten 1998; Levin 1999), 
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such that the fiducial accretion torque yV accl is balanced by r- 
mode angular momentum loss dJ/dt\ gr . 

The gravitational radiation reaction adds energy to the unsta- 
ble r-mode at a rate 
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is the canonical energy of the (I = 2,m = 2) r-mode (Fried- 
man & Schutz 1978; Owen et al. 1998). In a steady state 
all of this energy must be dissipated by viscous processes at 
a rate Wd = dE c /dt. In terms of the accretion luminosity, 
La = GMM/R = N accr flK, the dissipation rate is 
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The viscosity in the neutron star originates from several pos- 
sible sources. For normal npe matter, calculations of the vis- 
cous transport coefficients exist only at near-nuclear densities 
(Flowers & Itoh 1979). The components of such a core are 
strongly degenerate, and phase-space restrictions impart a char- 
acteristic T~ 2 dependence to the shear viscosity (Cutler & Lind- 
blom 1987). Compressing a fluid element of neutron star mat- 
ter causes it to emit neutrinos as the npe mixture reestablishes 
/3-equilibrium, so the bulk viscosity has an Urca-like T 6 depen- 
dence (Sawyer 1989). Another possibility for the viscosity is 
that it is caused by mutual friction in the neutron-proton su- 
perfiuid (Mendell 1991). In this case the viscous damping is 
independent of temperature. 

While the total amount of viscous dissipation Wd depends 
only on the assumption of a steady-state r-mode amplitude, the 
amount of heat actually deposited into the star depends on the 
nature of the damping. If the dominant viscous mechanism is 
bulk viscosity (i.e., for core temperatures T ;> 10 9 K), then the 
dissipated energy is released in the form of neutrinos, which 
promptly leave the star. The core temperatures of LMXBs are 
most likely less than 10 9 K, however, in which case the dissipa- 
tion mechanism is either shear viscosity or mutual friction. For 
both of these mechanisms, the heat Wd is deposited directly into 
the core of the star; we shall assume this to be the case in the 
rest of this paper. 

Levin (1999) first noted that, if the nucleons in the core are 
normal, the r-modes damped by shear viscosity are likely to be 
thermally unstable, at least for saturation amplitudes of order 
unity. The heating from the shearing motions decreases the vis- 
cosity, and so the r-mode amplitude increases, which heats the 
star even more. The result is a thermal and dynamical runaway. 
As envisaged by Levin (1999), the neutron star enters a limit 
cycle of slow spin-up to some critical frequency, at which the 
r-mode becomes unstable, followed by a rapid spin-down until 
the mode is once again damped. As the neutron star cools, ac- 
cretion again exerts a positive torque on the star, and the cycle 
repeats. Because the r-modes are present at a nonzero ampli- 
tude for only ~ 10~ 7 of the entire cycle's duration, it is unlikely 
that any of the known LMXBs harbor active r-modes and have 
normal fluid cores. 

For a superfluid core, where the damping is due to mutual 
friction (and hence independent of temperature), the neutron 
star can reach a state of three-fold equilibrium (Bildsten 1998; 
Levin 1999): the temperature is set by the balance of viscous 
heating and radiative or neutrino cooling, the r-mode's am- 
plitude is set by the balance of gravitational radiation back- 
reaction and viscous damping, and the spin is set by the balance 
of accretion torque and angular momentum loss to gravitational 
radiation. It is this scenario that we shall examine for existing 
evidence of r-mode spin regulation. 

While a neutron star accretes, its luminosity is dominated by 
the release of the infalling matter's gravitational potential en- 
ergy, La Rj 190MeV(M/f7ift), where mj is the average nucleon 
mass. Nuclear burning (either steady or via type I X-ray bursts) 
of the accreted hydrogen and helium generates an additional 
~ 5MeV per accreted nucleon. Most of this heat is promptly 
radiated away, however, and no more than a few percent dif- 
fuses inward to heat the interior (Fujimoto et al. 1984, 1987). 
Nuclear reactions in the deep crust (at p ;> 5 x 10 n gcm" 3 ) 
release about lMeV per accreted nucleon (Sato 1979; Blaes 
et al. 1990; Haensel & Zdunik 1990) and heat the crust directly 
(Brown & Bildsten 1998; Brown 2000). 

In addition to the crustal reactions, the viscous dissipation 



of r-modes constitutes another heat source in the neutron star's 
core. For a fiducial neutron star with M = 1 .4M Q and R = 10 km, 
equation (||) implies that W d /L A = 0.046(t/ sp i n /300Hz), or 
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This heating is very substantial, as it is much greater than the 
amount of nuclear heating from the crustal reactions. The 
prospects for detecting the effect of core r-mode heating in 
steadily accreting neutron stars are dim, unfortunately, as it is 
dwarfed by the accretion luminosity (which is a factor of 20 
brighter). The thermal emission from the neutron star is directly 
observable, however, if accretion periodically halts, as in the 
neutron star transients (the cooling timescale of the heated core 
is ~ 10 4 yr). While continued accretion at low levels between 
outbursts may contribute some of the quiescent luminosity (see 
Brown et al. 1998 for a discussion), the thermal emission from 
the hot crust of the neutron star is impossible to hide, and so 
observations of L q set an upper limit on the core temperature. 
The neutron stars in soft X-ray transients therefore offer the 
best prospects to look for evidence of viscous heating. In the 
next section we predict the quiescent luminosity L q that arises 
because of the r-mode heating and compare it to the observed 
luminosities of several neutron star transients. 

3. THE QUIESCENT LUMINOSITIES OF NEUTRON STAR 
TRANSIENTS 

The neutron star accretes fitfully, so the spin period and the 
r-mode amplitude oscillate about the equilibrium defined by the 
time-averaged accretion rate, (M) = t~ l J Mdt, where t r is the 
recurrence interval. Moreover, the timescale for viscous dissi- 
pation to heat the core is 
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where c p is the specific heat per baryon and (Wd) is the viscous 
heating averaged over an outburst/quiescent cycle. Because 
is much longer than the outburst recurrence time (typically of 
order years to decades), the core should remain fixed at the 
temperature set by the balance (over many outburst/quiescent 
cycles) between heating and cooling processes. We may there- 
fore compute the viscous dissipation using (M). Some simple 
estimates of the equilibrium core temperatures and the result- 
ing quiescent luminosities, for when both radiativ e an d neu- 
trino cooling are im portant, are presented first (§ 3.1). This 
is followed, in § 3.2, by detailed numerical calculations of the 
neutron star's thermal structure and a comparison (§ [Q| ) to ob- 
servations of several neutron star transients. 

3.1. Simple Estimates 

In a thermal steady state, the neutron star interior is cooled 
both by neutrinos emitted from the core and crust and by pho- 
tons emitted from the surface. To begin, we estimate the lu- 
minosity and the equilibrium core temperature set by balancing 
the heat deposited during an outburst/recurrence cycle, (Wd), 
with each cooling mechanism individually. First, if neutrino 
emission from the core is negligible (e.g., if the core is super- 
fluid and the Urea processes are exponentially suppressed), then 
all of the heat generated by viscous dissipation, (Wd), is con- 
ducted to the surface of the neutron star and escapes as thermal 
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radiation during quiescence. For the interior to be in a thermal 
steady state, the quiescent luminosity must then be 



L q w (W d ) =5.4 x 10 33 ergs _1 
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This estimate depends only on the assumption that neutrino 
emission is suppressed, and is independent of the crust micro- 
physics. 

As a check, we estimate the temperature of the neutron star 
core. In quiescence the atmosphere and crust come to resemble 
a cooling neutron star (Bildsten & Brown 1997; Brown et al. 
1998). For the temperature increase through the atmosphere 
and upper crust, we use the fit of Gudmundsson, Pethick, & 
Epstein (1983), 
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where Tb is the temperature at a fiducial boundary pb = 
10 10 gcm~ 3 . Equating L 7 with L q from equation (||) gives an 
estimate of the temperature in the upper crust, 
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To relate Tb to the core temperature T c , we use approximate ana- 
lytic expressions (Brown 2000; eqs. [22] and [23]) for the crust 
temperature to obtain 
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where we have neglected the luminosity due to crustal nuclear 
reactions. Substituting from equation ( |ll| ) for Tb, we obtain the 
core temperature in the absence of neutrino emission, 
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where the scalings for (M) and ^ sp i n are obtained by dropping 
the second term on the right in equation (12). This estimate 
agrees quite well with the detailed calculations described in 



§ 3.2 



The core neutrino emissivity is, for modified Urea processes 
(Shapiro & Teukolsky 1983), L^ lca w 7.4 x 10 31 (r c /10 8 K) 8 , 
multiplied by a superfluid reduction factor that goes roughly 
as exp(-A/kT c ), where A is the superfluid gap energy 
(Yakovlev & Levenfish 1995). For A > kT c the net Urea neu- 
trino luminosity is much less than L q , so that equation (Q) 
is self-consistent. Neutrino emission from crust neutrino 
bremsstrahlung (Kaminker et al. 1999) at the temperature 7j 
(eq. [pl]|) is also not significant, although at higher (M) it is 
competitive with radiative cooling. Hence for accretion rates 
typical of neutron star transients, the majority of the deposited 
heat is conducted to the surface, and equation (^) provides a 
robust estimate of the radiative luminosity of the star. 

Alternatively, if core neutrino emission is not suppressed 
(i.e., the nucleons are not superfluid), then modified Urea pro- 
cesses are the dominant coolant and L^ rca rs Wd- In this case the 
core temperature is 
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and is smaller than if the core were superfluid. A colder core 
implies a dimmer thermal luminosity from the surface. In order 



to estimate L q , we write Wd 



, where L q = L 7 (7},), 



and Tb is related to T c by equation (|12|). Under the assump- 
tion that L q <C L^ rca , the solution of the resulting transcendental 
equation is 
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where we obtain the scalings for (M) and f sp ; n by dropping the 
second term on the right-hand side of equation (|l2|). In this 
case L q is less than W d and L^ rca , so our assumption that core 
neutrino emission is the dominant coolant is self-consistent. 

These estimates neglect cooling from other neutrino- 
producing mechanisms, such as neutrino bremsstrahlung in the 
crust. Moreover, the core neutrino emissivity depends on the 
local proper temperature, which increases towards the center of 
the star because of the gravitational redshift. We now describe 
our detailed calculations, which take these effects into account. 

3.2. Numerical Calculations 

To calculate the expected quiescent luminosities of accret- 
ing neutron stars, we compute hydrostatic neutron star mod- 
els by integrating the post-Newtonian stellar structure equations 
(Thorne 1977) for the radius r, gravitational mass m, potential, 
and pressure with the equation of state AV18+<5v+UIX* (Ak- 
mal, Pandharipande, & Ravenhall 1998), as described in Brown 
(2000). With the hydrostatic structure specified, the luminosity 
L and temperature T are found by solving the entropy and flux 
equations (Thorne 1977), 
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Here e r and are the nuclear heating and neutrino emissiv- 
ity per baryon, n is the baryon density, and K is the thermal 
conductivity. The potential $ appears in the time-time compo- 
nent of the metric as e*/ e (it governs the redshift of photons 
and neutrinos; Misner, Thorne, & Wheeler 1973). We neglect 
in equation ( fl~^ ) terms arising from compressional heating, as 
they are of order TAs(M/M) (Fujimoto & Sugimoto 1982), s 
being the specific entropy, and are negligible throughout the de- 
generate crust and core (Brown & Bildsten 1998). We do not 
include heating from nuclear reactions in the deep crust. This 
has the effect of underestimating slightly (by <J 10%) the qui- 
escent luminosity of the neutron star transient. Equations ( Jig ) 
and ( |l7| ) are integrated outwards to a density pb = 10 10 gcrrf^\ 
We there impose a boundary condition relating L and T with the 
fitting formula of Potekhin, Chabrier, & Yakovlev (1997) for a 
partially accreted crust. By incorporating a parameter describ- 
ing the depth of a light element (H and He) layer, this formula 
differs from that of Gudmundsson et al. (1983), which we used 
for our simple estimates (§ 3.1). We set the depth of this light 



element layer to where the density is w 10 gem , which is 
roughly where the accreted material burns to heavier elements 
(Hanawa & Fujimoto 1986). 
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The high thermal conductivity of the neutron star's core in- 
sures that it is very nearly isothermal, regardless of the detailed 
dependence of the heating rate e r on the radius. For the core 
temperatures typical of LMXBs, bulk viscosity is unimportant, 
so we assume that the heating is from ordinary shear viscos- 
ity. The rate per unit volume is just 2rjSa ah Sa* b , where ?y is the 
shear viscosity and 5a a b is the kinematic shear. If we neglect 
the dependence of shear viscosity on density (since the density 
is approximately constant in the neutron star's core), this rate 
is just proportional to r 2 for an (/ = 2, m = 2) r-mode (Lindblom 
et al. 1998). Hence we take e r oc r 2 , and normalize it so that 
the heating rate, when integrated over the core, satisfies equa- 
tion (g). 

The microphysics used to integrate equations ( |l6j ) and ( |l7| ) 
is fully described in Brown (2000), so here we just highlight 
two modeling uncertainties. First, standard calculations pre- 
sume that the neutron star's crust is a pure lattice, and hence the 
conductivity is dominated by electron-phonon scattering. Over 
the lifetime of an LMXB, however, the neutron star can eas- 
ily accrete enough matter to replace its entire crust (requiring 
about 0.01M Q ). The accreted crust is formed from the products 
of hydrogen and helium burning and is likely to be very impure 
(Schatz et al. 1999). A lower conductivity from impurities low- 
ers the surface temperature, and hence the quiescent luminosity, 
for a given core temperature. We model the low thermal con- 
ductivity of a very impure crust by using electron-ion scattering 
(Haensel, Kaminker, & Yakovlev 1996) throughout the crust. 

The second modeling uncertainty is the superfluid transi- 
tion temperatures, for which estimates vary widely (see Tsu- 
ruta 1998, and references therein). When the core temperature 
is much less than the superfluid transition temperature, emissiv- 
ity from Cooper pairing is unimportant (Yakovlev, Kaminker, 
& Levenfish 1999) and the superfluidity suppresses the neu- 
trino emission by roughly the Boltzmann factor exp(-A/lcBT), 
for a superfluid gap energy A. We perform our calculations for 
two models, one with superfluidity parameterized as in Brown 
(2000), with a typical gap energy A ~ 0.5 MeV, and another 
model with a normal core, A = 0. 

Figure 1 demonstrates the thermal structure of such a neu- 
tron star with a time-averaged accretion rate (M) = 2.4 x 
10 -1 ' Mq yr - ' (the rate inferred for Aql X-l) for a spin fre- 
quency of 275 Hz (solid lines) and 549 Hz (dotted lines). If the 
core is superfluid (the upper pair of curves), then the neutrino 
luminosity from crust bremsstrahlung (region leftward of the 
vertical dot-dashed line) is roughly comparable to the photon 
luminosity. In contrast, if the core were normal (so that the 
modified Urea processes were unsuppressed) but the viscosity 
remained independent of temperature (so that a thermal steady 
state could be reached), then only about 10% of the heat gener- 
ated in the core would be conducted to the surface. The rest of 
the heat is balanced by modified Urea neutrino emission. The 
core temperature and fraction of viscous hea t co nducted to the 
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FIG . 1 — The thermal structure of a neutron star accreting at a time-averaged 
rate of 2.4 X 10 -11 Mq yr~' (e.g., Aql X-l), for two different spin frequencies: 
275 Hz (solid lines) and 549 Hz (dotted lines). The upper pair of curves are for 
a superfluid core (region rightward of the vertical dot-dashed line); the lower 
pair, for a normal core. 



3.3. Comparison to Observed Transients 

A superfluid core is cooled mainly by conduction of heat 
to the surface, at least until the interior temperature is high 
enough to activate crust neutrino bremsstrahlung. Figure 2 
shows the expected quiescent luminosity L q as a function of 
(M) for this case, with a range (shaded region) of rotation fre- 
quencies 200 Hz < f < 600 Hz. The inferred (M) and L q for 
several neutron star transients are also plotted (squares) for 
comparison. With the exception of EXO 0748-676 2 , the neu- 
tron star transients with measured quiescent luminosities are 
too dim, by a factor of 5-10, to be consistent with viscous heat- 
ing of the magnitude assumed here. We must conclude, then, 
that either the accretion torque is much less than (M)(GMR) 1 ' 2 , 
or that a steady-state r-mode does not set their spin. 

The quiescent luminosities for Aql X-l, Cen X-4, and 
4U 1608-522 use the bolometric corrections appropriate for a H 
atmosphere spectrum (Rutledge et al. 1999a); L q for the Rapid 
Burster is from Asai et al. (1996a). We infer the time-averaged 
accretion rate from (M) w (t„/t r )(L /GMR~ l ), where 

2 EXO 0748-676 is likely to accrete during quiescence, as suggested by observations (with ASCA) of variability on timescales > 1000s (Corbet et al. 1994; Thomas 
et al. 1997). 
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FIG 2. — Quiescent luminosities as a function of time-averaged accretion rate 
(M). The heating from viscous dissipation of the r-mode is from eq. (^), and 
the neutrino emission from the core is suppressed by nucleon superfluidity. 
The shaded region corresponds to rotation frequencies between 200 Hz (lower 
curve) and 600 Hz (upper curve). Neutrino cooling from crust bremsstrahlung 
is important rightward (i.e., at higher (M)) of the knee in the shaded region. 
Also shown are the inferred quiescent luminosities and time-averaged accre- 
tion rates for several neutron star transients. 

t a and L a are the outburst duration and luminosity and the dis- 
tances are taken from Chen, Shrader, & Livio (1997). 3 Out- 
burst fluences for Aql X-l and the Rapid Burster are accurately 
known (RXTEI AW-Sky Monitor public data); for the remaining 
sources (M) is estimated from peak luminosities and outburst 
rise and decay timescales (Chen et al. 1997). 

Our estimates for (M) depend on the inferred source dis- 
tance. When most of the r-mode heating Wd is conducted to 
the surface, however, as in the superfiuid core case for (M) <^ 
10 _1 1 Mq yr -1 , the predicted quiescent luminosity is L q m Wd oc 
(M) (see eqs. [^] and [[7]]), and hence depends on the source 
distance in the same way as does (M). Therefore, our com- 
parison of L q predicted from r-mode heating and the quiescent 
luminosity actually observed is independent of distance. In this 
regime, the relation between L q and (M) is also independent of 
the microphysics in the crust. 

As shown by Levin (1999), the temperature dependence of 
viscosity in a normal fluid likely prevents a steady-state r-mode. 
For comparison, however, we plot in Figure 3 the case where 
modified Urea neutrino emission from the core is allowed (as it 
would be in a normal fluid) but the r-mode amplitude is steady, 
i.e., we assume that a thermogravitational runaway has some- 
how been avoided. As a result, neutrino emission efficiently 
cools the core, and so the radiative luminosity L q is less for 
a given (M). For this case, with the exception of Cen X-4, 
the neutron star transients have quiescent luminosities roughly 
consistent with that predicted. Because there is a characteristic 
core temperature, namely, that at which neutrino cooling equals 
radiative cooling, the relation between (M) and L q is no 
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FIG. 3 — The same as Fig. 2, but for a normal core. We also show (thin 
dotted lines) L q ((M)) for a neutron star with a crust of light elements at 
p < 10 I0 gcnr 3 . 

longer independent of distance, unlike the case shown in Fig- 
ure 2. The knee in the shaded region is where the neutrino and 
photon luminosities are comparable. Rightward of this knee 
{(M) J> 10~ 12 MQyr _1 ) neutrino cooling prevents the core tem- 
perature, and hence the photon luminosity, from rising rapidly 
with increasing (M). Should the crust have a higher conduc- 
tivity (e.g., if it were more pure) than we have assumed here, 
then the shaded region rightward of the knee would move up- 
wards, i.e., the predicted L q would be even higher. To illustrate 
this we computed L q ((M)) using the L 7 (7i) relation for a crust 
composed of light elements (and having a higher conductivity) 
for densities less than p\, = 10 10 gcm -3 (dotted lines). 

It should be noted that the actual thermal radiation from a 
neutron star's surface is in general less than the observed quies- 
cent luminosity, since other emission mechanisms are possible, 
such as accretion via a low-efficiency advective flow (Narayan, 
McClintock, & Yi 1996) or magnetospheric emission (Cam- 
pana et al. 1998a). Evidence for other, non-thermal emission 
processes are the hard power-law tails observed from Cen X- 
4 (ASCA; Asai et al. 1996b) and Aql X-l (BeppoSAX; Cam- 
pana et al. 1998b). In addition, variability on timescales of a 
few days has been observed from Cen X-4 (van Paradijs et al. 
1987; Campana et al. 1997). As a result, a plot showing ther- 
mal emission (as opposed to observed L q ) would have the data 
points shifted downward in Figures 2 and 3. In other words, 
the quiescent luminosity inferred from observations is likely to 
overestimate the actual thermal emission from the neutron star. 
This strengthens our conclusion regarding the incompatibility 
of steady-state r-mode heating with the observations. 

There are stronger neutrino emission mechanisms possible 
than modified Urea and crust bremsstrahlung. Recently, there 
has been renewed interest in the direct Urea process (Lattimer 
et al. 1991), which is allowed if the proton fraction exceeds 
0. 148 or if hyperons are present (Prakash et al. 1992). Other ex- 
otic mechanisms may be possible, including pion condensates 



3 Recent observations (Callanan, Filippenko, & Garcia 1999) resolved the optical counterpart of Aql X-l into two objects. We use the distance estimate (2.5kpc) 
of Chevalier et al. (1999), which accounts for the interloper star. 
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(Umeda et al. 1994), kaon condensates (Brown et al. 1988), or 
quark matter (Iwamoto 1982). The exotic mechanisms have 
the same temperature dependence as the direct Urea (oc T 6 ) 
but are weaker. Should any of these enhanced processes oc- 
cur, the core will be much colder, and the heat radiated from 
the surface much weaker, than in the calculations here. For 
example, balancing the viscous heating with neutrino emis- 
sion from apion condensate, ss 2.0 x 10 39 (r/10 8 K) 6 ergs _1 
(Shapiro & Teukolsky 1983), implies that T c k, 1.2 x 
lO 7 ((M)/lO- 11 M yr 1 ) 1 / 6 K, and, from equations @ and 
(H), that L q « 6.0 x lO 3O ((M}/lO- 11 M yr 1 ) a45 ergs- 1 . This 
is much dimmer than that observed. Of course, it is possible 
that superfluidity reduces L v such that the core temperature is 
just enough to explain the observed quiescent emission. It is 
difficult, however, to arrange all of the sources to obey such a 
relation. 

4. CONCLUSIONS 

Using the assumption that the accretion torque is balanced 
by angular momentum loss from gravitational radiation by an r- 
mode pulsation of constant amplitude, we find that the expected 
quiescent luminosities of the neutron star X-ray transients, for 
rotation rates of 200-600 Hz, are characteristically brighter than 
those observed. Reconciling the observations with the presence 
of r-mode heating requires that neutrino emission from the core 
be unsuppressed, as for a normal core. In this case, however, 
the r-mode is thermally unstable and cannot remain at a con- 
stant amplitude, unless some mechanism prevents a runaway. 
It therefore seems unlikely that the spin frequency of Aql X- 
1 is a signature of a steady-state core r-mode pulsation. We 
note, however, that the same conclusion cannot be drawn for 
the bright, persistent LMXBs (such as Sco X-l, which could be 
detected by gravitational wave experiments soon to be opera- 
tional). Uncertainties in the nuclear burning and the accretion 
luminosity cannot constrain the surface thermal luminosity to 
within <; 5%, which is necessary to differentiate the r-mode 
heating from the accretion luminosity. 

In addition to Aql X-l, there is one other neutron star tran- 
sient which is known to be spinning rapidly, and that is the 
401 Hz accreting pulsar (Wijnands & van der Klis 1998) in 
the transient SAX J1808.4-3658 (in't Zand et al. 1998). This 
source has not yet been detected in quiescence. Given a recur- 
rence interval of 1 .5 yr, an outburst duration of w 20 day, and an 
outburst accretion rate of w 3 x 10~ 10 MQyr"' (in't Zand et al. 
1998), we expect a quiescent luminosity > 5.8 x 10 33 ergs _1 if 
the core is superfluid and an active r-mode pulsation balances 
the accretion torque in this system. This L q corresponds to 
an unabsorbed flux (4kpc distance) of 3 x 10~ 12 ergcm" 2 s" 1 , 
which is about ten times the flux expected if the only heat 
source were crust nuclear reactions (Brown et al. 1998). Fu- 
ture ASCA, Chandra, and XMM observations will assist in con- 
straining the viscous damping present. The luminosity from 
the viscous damping is much larger than the expected magneto- 
spheric emission (Becker & Triimper 1997), and so interpreta- 
tion of the spectrum should be unambiguous in the absence of 
accretion onto the neutron star's surface. 

At (M) <; 10 _11 Mq yr" 1 , all of the viscous heating in the core 
is radiated from th e ne utron star's surface during quiescence. 
As noted in section 3.3, the relation between L q ((M)) then de- 
pends only on the accretion torque, and not on the source dis- 
tance and crust microphysics. Chandra and XMM are ideally 



suited for a study of a population of low-luminosity neutron 
stars, which offer excellent prospects for a clean determination 
of the amount of viscous heating present. 

At higher (M), for which neutrino cooling from the crust 
contributes to balancing the viscous heating, the quiescent lu- 
minosity depends on the crust microphysics. In our calculations 
we assume that the neutron crust is very impure and hence used 
thermal conductivity dominated by electron-ion collisions. If 
the conductivity of the crust is higher than we have assumed 
(e.g., if the crust is a pure lattice), then the predicted quies- 
cent luminosity L q would be even higher than that plotted in 
Figures 2 and 3. In addition, we underestimated the predicted 
L q by neglecting the effect of direct heating of the neutron star 
crust by nuclear reactions occurring near neutron drip (Brown 
et al. 1998). Moreover, taking into account the possibility that 
non-thermal emission contributes to the observed quiescent lu- 
minosity further widens the gap between the observed L q and 
that inferred from the r-mode spin regulation hypothesis. All of 
these effects further strengthen our conclusions. 

If the r-mode is not in steady state, then there remain several 
possibilities: either the superfluid viscosity is so strong that it 
suppresses the r-mode instability entirely, or the mode satura- 
tion amplitude is so small that it is unimportant at all the spin 
frequencies observed, or else the neutron star is in a limit cycle 
(Levin 1999) of spin-up to some critical frequency, followed 
by rapid spin-down and heating. A detailed study of the spin 
evolution is necessary to determine if the spin periods of the 
neutron stars are consistent with such a scenario. In particular, 
it remains an open question as to whether one should expect to 
observe a population of slowly spinning neutron stars with low- 
mass companions, such as Her X-l, 4U 1626-67, and GX 1+4. 

The study of r-modes in neutron stars is rapidly evolving in 
response to the interest aroused in the general relativity com- 
munity. While this paper was being refereed, several theoreti- 
cal developments occured that are relevant for this study. First, 
Lindblom & Mendell (1999) showed that unless the superfluid 
entrainment parameter assumes a very special value, superfluid 
mutual friction is not competitive with gravitational radiation 
for the r-mode amplitude evolution. There is therefore a con- 
flict between theory and experiment: while theoretical calcula- 
tions show that r-modes in superfluid neutron stars should be 
excited, the observations discussed in this paper are direct evi- 
dence against the r-modes having a sufficient steady amplitude 
to limit the spin of the neutron star, and the clustering of LMXB 
spin frequencies argues against an recurrent instability. This 
contradiction is likely resolved by consideration of the presence 
of a solid crust (Bildsten & Ushomirsky 1999), which dramat- 
ically enhances the dissipation rate and damps the r-modes for 
typical core temperatures and spin frequencies of LMXBs. The 
findings presented in this paper lend observational support to 
that conclusion. 
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supported by NASA via grant NAGW-45 17. E.FB is supported 
by a NASA GSRP Graduate Fellowship under grant NGT5- 
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and John Hertz Foundation. 



8 



Brown and Ushomirsky 



REFERENCES 



Akmal, A., Pandharipande, V. R., & Ravenhall, D. G. 1998, Phys. Rev. C, 58, 
1804 

Andersson, N. 1998, ApJ, 502, 708 

Andersson, N., Kokkotas, K. D., & Schutz, B. F. 1999a, ApJ, 510, 846 
Andersson, N., Kokkotas, K. D., & Stergioulas, N. 1999b, ApJ, 516, 307 
Asai, K., Dotani, T., Kunieda, H., & Kawai, N. 1996a, PASJ, 48, L27 
Asai, K., Dotani, T., Mitsuda, K., Hoshi, R., Vaughan, B., Tanaka, Y., & Inoue, 

H. 1996b, PASJ, 48, 257 
Becker, W. & Triimper, J. 1997, A&A, 326, 682 

Bhattacharya, D. 1995, in X-ray Binaries, ed. W. H. G. Lewin, J. van Paradijs, 

& R. E. Taam (Cambridge University Press), 233-251 
Bildsten, L. 1998, ApJ, 501, L89 
Bildsten, L. & Brown, E. F. 1997, ApJ, 477, 897 

Bildsten, L. & UshamilSky, G 19Q9 ApT in pregg prpprinf avnilahle- gcfrr.- 

ph/99 1 1 1 55 

Blaes, O., Blandford, R., Madau, P., & Koonin, S. 1990, ApJ, 363, 612 



nn T IQQQ Phv 



n 



Prpprmf nvnilnhlp ^ gf- 



Brady, P. & Cp 

qc/9812014 1 - 

Brown, E. F. 2000, ApJ, in press. Preprint available: astro-ph/9910215 

Brown, E. F. & Bildsten, L. 1998, ApJ, 496, 915 

Brown, E. F, Bildsten, L., & Rutledge, R. E. 1998, ApJ, 504, L95 
Brown, G, Kubodera, K., Page, D., & Pizzochero, P. 1988, Phys. Rev. D, 37, 
2042 

Callanan, P. J., Filippenko, A. V., & Garcia, M. R. 1999, IAU Circ, 7086 
Campana, S., Colpi, M., Mereghetti, S., Stella, L., & Tavani, M. 1998a, 
A&A Rev., 8, 279 

Campana, S., Mereghetti, S., Stella, L., & Colpi, M. 1997, A&A, 324, 941 
Campana, S., Stella, L., Mereghetti, S., Colpi, M., Tavani, M., Ricci, D., Dal 

Fiume, D., & Belloni, T. 1998b, ApJ, 499, L65 
Chen, W., Shrader, C. R., & Livio, M. 1997, ApJ, 491, 312 
Chevalier, C, Ilovaisky, S. A., Leisy, P., & Patat, F. 1999, A&A, 347, L51 
Corbet, R. H. D., Asai, K., Dotani, T., & Nagase, F. 1994, ApJ, 436, L15 
Cutler, C. & Lindblom, L. 1987, ApJ, 314, 234 
Flowers, E. & Itoh, N. 1979, ApJ, 230, 847 
Friedman, J. L. & Morsink, S. M. 1998, ApJ, 502, 714 
Friedman, J. L. & Schutz, B. F. 1978, ApJ, 222, 281 

Fujimoto, M. Y., Hanawa, T., Iben, I., & Richardson, M. B. 1984, ApJ, 278, 
813 

— . 1987, ApJ, 315, 198 

Fujimoto, M. Y. & Sugimoto, D. 1982, ApJ, 257, 291 
Ghosh, P. & Lamb, F. K. 1979, ApJ, 234, 296 

Gudmundsson, E. H., Pethick, C. J., & Epstein, R. I. 1983, ApJ, 272, 286 
Haensel, P., Kaminker, A. D., & Yakovlev, D. G. 1996, A&A, 314, 328 
Haensel, P. & Zdunik, J. L. 1990, A&A, 227, 431 
Hanawa, T. & Fujimoto, M. Y. 1986, PASJ, 38, 13 

in't Zand, J. J. M., Heise, J., Muller, J. M., Bazzano, A., Cocchi, M., Natalucci, 
L., & Ubertini, P. 1998, A&A, 331, L25 



Iwamoto, N. 1982, Ann. Phys., 141, 1 

Kaminker, A. D., Pethick, C. J., Potekhin, A. Y, Thorsson, V., & Yakovlev, 

D. G. 1999, A&A, 343, 1009 
Lattimer, J. M., Pethick, C. J., Prakash, M., & Haensel, P. 1991, 

Phys. Rev. Lett., 66, 2701 
Levin, Y. 1999, ApJ, 517, 328 

Linrlhlnm , T & Mendell, G. 1999, Phys. Rev. D, submitted, preprint available: 



M. 1998, Phys. Rev. Lett., 80, 4843 



gr-qc/9909084 
Linawom, L., uwen, B. J., & Morsink, S. 
Mendell, G. 1991, ApJ, 380, 530 
Miller, M. C. 1999, ApJ, 515, L77 

Misner, C. W., Thorne, K. S., & Wheeler, J. A. 1973, Gravitation (New York: 

W. H. Freeman and Company) 
Narayan, R., McClintock, J. E., & Yi, I. 1996, ApJ, 457, 821 
Owen, B. J., Lindblom, L., Cutler, C, Schutz, B. F, Vecchio, A., & Andersson, 

N. 1998, Phys. Rev. D, 58, 084020 
Papaloizou, J. & Pringle, J. E. 1978, MNRAS, 184, 501 
Potekhin, A. Y, Chabrier, G., & Yakovlev, D. G. 1997, A&A, 323, 415 
Prakash, M., Prakash, M., Lattimer, J. M., & Pethick, C. J. 1992, ApJ, 390, L77 
Rutledge, R. E., Bildsten, L., Brown, E. F—EaalflK, G G A Zavlin, V. E. 1999a, 

ApJ, 529, in press, preprint available: istro-ph/990931S 

— . 1999b, ApJ, 5 14, 945 1 

Sato, K. 1979, Prog. Theor. Physics, 62, 957 
Sawyer, R. F. 1989, Phys. Rev. D, 39, 3804 

Schatz, H, Bildsten, L., Cumming, A., & Wiescher, M. 1999, ApJ, 524, 1014 
Shapiro, S. L. & Teukolsky, S. A. 1983, Black holes, white dwarfs, and neutron 

stars (New York: John Wiley & Sons, Inc.) 
Strohmayer, T. E., Zhang, W., Swank, J. H, Smale, A., Titarchuk, L., Day, C, 

& Lee, U. 1996, ApJ, 469, L9 
Thomas, B., Corbet, R., Smale, A. P., Asai, K., & Dotani, T. 1997, ApJ, 480, 

L21 

Thorne, K. S. 1977, ApJ, 212, 825 
Tsuruta, S. 1998, Phys. Rep., 292, 1 

Umeda, H, Nomoto, K, Tsuruta, S., Muto, T, & Tatsumi, T. 1994, ApJ, 431, 
309 

van der Klis, M. ] 999, in Proceedings of the Third William Fairbank Meeting, 

available from astro-ph/98 1 2395 
van Paradijs, J., Verbunt, K, Shater, R. A., & Arnaud, K. A. 1987, A&A, 182, 

47 

van Paradijs, J. & White, N. 1995, ApJ, 447, L33 
Wagoner, R. V. 1984, ApJ, 278, 345 
White, N. E. & Zhang, W. 1997, ApJ, 490, L87 
Wijnands, R. & van der Klis, M. 1998, Nature, 394, 344 
Yakovlev, D. G, Kaminker, A. D., & Levenfish, K. P. 1999, A&A, 343, 650 
Yakovlev, D. G. & Levenfish, K. P. 1995, A&A, 297, 717 
Zhang, W., Jahoda, K, Kelley, R. L., Strohmayer, T. E., Swank, J. H, & Zhang, 
S. N. 1998, ApJ, 495, L9 



